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ON THE GAS DYNAMICS OF A ROTATING IMPELLER* 
By A. Busemann 


INTRODUCTION 


Centrifugal pumps and turbines may sometimes be treated by using 
plane incompressible flow where the fluid springs from a source or runs 
into a sink in the center of the impeller. However, application of 
conformal - mapping is not so simple for this type of flow as for recti- 
linear flow, because a steady flow about the blade sections exists only 
in a coordinate system rotating with the blade. In a coordinate 
system at rest, the velocity field outside of the singularities is free 
from sources and rotation, but the impeller rotates in it. By applying 
appropriate distributions of sources, sinks, and vortices inside the 
blade contour, a flow can be found which does not pass through the 
rotating blade contour. The action of the flow on these singularities 
inside the blades creates the torque. 

For blowers with high circumferential velocities, the same consider- 
ations would have to be extended by the requirements of a compressible 
flow. Though there exist corresponding laws of forces on sources, sinks, 
and vortices, the conditions become considerably more involved because 
the freedom from vortices in the regular domain of the flow concerns the 
velocity field whereas freedom from sources exists only In steady flow and 
involves the stream-density field (the product of velocity and gas density) . 
It becomes impossible to superimpose two fields free ipcan sources and 
vortices into a new field which is free from sources and vortices; 
furthermore, in a compressible flow one can operate with distributed 
singularities only, since a certain region around point sources and point 
vortices Is void of any velocity field. For these reasons, the laws of 
forces on singularities have only a very limited range of application in 
gas dynamics . 

It is known that the velocity fields of incompressible flow and 
magnetic fields are similar with respect to the distribution of the 
vector and to the field energy. Accordingly, one should anticipate the 
forces on corresponding singularities to be equal so that a displacement 
of the singularities produces the amount of work required for the change 
in field energy connected with it. Contrary to this expectation the 
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hydrodynamic and the magnetic forces are always opposite in sign. The 
result is that, on the one hand, the forces on magnetic sources and sinks 
(north and south poles) and. those on hydrodynamic vortices are such that 
their displacement work balances the field energy. On the other hand, 
for the displacement of hydrodynamic sources and sinks and of magnetic 
vortices (around conductors of electric currents), the opposite sign of 
the displacement work would be required to balance the field energy. 

The difference is restored for the magnetic vortex according to Faraday's 
induction law by the well-known fact that the electric current flowing 
in- the conductor is opposed, during any displacement, by induced elec- 
tromotive forces, the overcoming of which requires a supply of electrical 
energy. Likewise, additional energy must be supplied during displace- 
ments of the hydrodynamic sources and sinks} in this case the energy 
supply is caused by higher pressures at the locatipn of the sources and 
lower pressures at the location of the s ink s compared to the pressures 
of the steady flow. These pressure differences are precisely what 
constitutes the pressure heads or pressure differentials for rotating 
machinery. 

For plane flow, one may regard instead of the field of the stream- 
lines the field of the potential lines which are orthogonal to the 
streamlines . In the case of incompressible flow, this field is in the 
regular domain free from sources and vortices; however, sources nnri 
vortices are interchanged, since a streamline source represents a 
potential- line vortex and vice versa. - Comparing plane potential lines 
with the plane field of magnetic lines of force, one finds now perfect 
agreement including the sign of the forces and the induction law. If 
one wants to determine' the additional pressure difference between two 
points of the field which alters Bernoulli's pressure difference of the 
steady flow, one must draw a line connecting these points and observe 
on it the variation with time In the number of the potential lines which 
intersect this connecting line. The connecting line must not pass over 
source points of the potential lines . This is the same rule which 
applies for the determination of induced electromotive forces . In order 
to ascertain whether the connecting line has moved over a source, these 
sources cannot be allowed to appear and disappear. For magnetic sources, 
this prerequisite is obviously satisfied by the fact that a north pole 
can be created only by separating a north and a south pole of equal 
strength. The sources of the potential lines represent- vortices of the 
streamlines and it follows from Helmholtz ' vortex theorems that, in two- 
dimensional flow, a vortex rotating clockwise can be produced only by 
separating a vortex rotating clockwise from another one rotating counter- 
clockwise. For these plane fields, the similarity is therefore very far 
reaching. For magnetic and hydrodynamic- fields in space, the similarity 
is subject to limitations because lines of forces remain lines in space 
whereas the potential lines change to potential surfaces insofar as 
unique surfaces orthogonal to the streamlines exist. Generally, one is 
therefore limited to correlate only the 'hydrodynamic velocity field and 
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the magnetic force field. If one identifies velocity and magnetic-field 
intensity, then the source matches the magnetic pole and the hydrodynamic 
vortex matches the wire through which current flows; with respect to the- 
induction law, however, it is just the wire with electric current and the 
hydrodynamic source which can he compared. 

Whereas the laws about forces on sources, sinks, and vortices almost 
maintain their form in transfer from hydrodynamics to gas dynamics hut 
lose their major field of application, the hydrodynamical induction law 
must correspond to a s imila rly formulated gas -dynamical induction law 
n.nri at the same time maintain to some extent its applicability, since 
it deals only with the pressure difference of unsteady flow as compared 
to steady flow. It is just this pressure difference which is related to 
the pressure rise in blowers . Thus , the author will give below the 
general derivation of the gas -dynamical induction law as he presented it 
for the first t ime in the s imme r of 1936 in a colloquium on gas dynamics 
at the DVL under the chairmanship of E. Trefftz. In a second part, he 
will show that the torque at the impeller in nonviscous compressible 
flow, even for circumferential velocities below sonic velocity, does not 
result from the effective pressure rise alone as in nonviscous incom- 
pressible flow, but that, on the contrary, even without effective pressure 
rise, energy quantities may be radiated from the rotating impeller. 


I. TEE PRESSURE RISE OF A CENTRIFUGAL BLOWER 
1 . Derivation of the Gas-Dynamical Induction Law 


In the revision of the section "Hydrodynamics " for the 8 th edition 
of A. Foppl, "Vorlesungen uber technische Mechanik" (Lectures on technical 
mechanics) , volume IV, page 417 , I derived the "hydrody n a mi cal induction 
law." If one limits the application to a loss-free gas flow of constant 
entropy, the derivation may be transferred directly to gases. 


For a nonviscous gas, free from gravity, of pressure p, density p, 
and velocity W with the components u, v, w in the directions of 
the spatial coordinates x, y, z, one obtains, in dependence on the 
t ime t, the following equations of motion: 


dp _ du dp _ dv dp _ dw 

dx ^ dt* dy ^ dt* dz dt 


( 1 ) 


For constant entropy s there applies for the enthalpy i of the gas: 



( 2 ) 
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By substitution of this relationship into equation (l) there result the 
following equations: 


_ 8i = du _ „ dy _ = d£ (3) 

dx dt dy dt dz dt 

For a certain time t = t Q one can combine these three partial deriva- 
tives into the following total differential of the enthalpy in space: 

-di = — dx+ — dy + — dz ( 4 ) 

dt dt dt 

If one introduces the convention that the components dx, dy, dz of a 
line element in space (which are used in equation ( 4 ) only for the 
time t = t D ) shall be for all times, the components of. a line element 
attached to the gas which connects the adjacent points and G2 

moved with the gas, the following transformation of equation ( 4 ) holds: 

-di = — (udx+vdy + w dz) - u du - v dv - w dw (5) 

dt 


since in this case the interchange of the differentiation 

— (dx) = d/— ^\= du, etc., is valid. If one places the points Gj_ and G2 
dt \dt/ 

attached to the gas farther apart, one may integrate the differentials 
indicated in equation (5) along a line from G]_ to G2 moved with the 

gas, so that the following difference in enthalpy is then obtained: 



If one introduces, instead of the points Gj_ and G2 attached to the 
gas, the pointB fixed in space P 1 and P 2 , which at the time t = t 0 
coincide with the former, one can prolong the line between G]_ and G2 
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with those gas points sweeping over the points P}_ and P2 after the 
time t = t Q (fig. l) . The integration along this line, which is like- 
wise attached to the gas (and is therefore determined in the development 
with time) hut connects the fixed points P-j_ anfl P 2 , includes beyond 

the integral required in equation (6) contributions which result from the 
shift of the integration limits and have the value dt. If 

one subtracts the latter, one obtains the following relation: 



(u dx + v dy + w dz) 


- i (uj2 + vj_2 + w-j_2) + 


| (u g 2 + v 2 2 + w £ 2) 


or 


*2 + \ (^2 2 + v 2 2 + w 2 2 ^ “ [^1 + 2 ^ U l 2 + v l 2 + V 1 2 |J = 

r* 2 (u 

Jp i 


_d_ 

dt 


dx + v dy + w dz) 


(7) 


The right side of the equation represents the induced pressure rise which 
is generated between the points P]_ and P 2 . Note that one deals here 

not with a partial differentiation with respect to time but with a total 
differentiation along a line moved with the gas. The selection of the 
line at the time ‘t Q is arbitrary. If this line, once selected, is not 
moved exactly with the gas, there originates an error which is proportional 
to the vortices located between the line moved exactly with the gas and the 
wrong connecting line. 


2 . Applications of the Induction Law 

In a steady gas flow free from vortices, the right side of equation (7) 
disappears because in the first place the integral, due to the freedom 
from vorticity, is independent of the path and equals the potential differ- 
ence $2 - and in the second place, this potential difference, due to 

the steady state, is independent of the time. 
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For a steady gas flow with vortices, the right side of equation (7) 
disappears only for points Pp and P2 which lie on the same stream- 
line. If one draws the connecting line from Pp to P2 along this 
streamline, all its points remain on this streamline later on, too. The 
entire integration path varies with time only in that it extends beyond 
the downstream point P2 in a flat loop which makes, however, no contri- 
bution to the integral. The value of the integral on the remaining piece 
from Pp to P2 is again Independent of the time so that the right side 
of equation (7) disappears. The vanishing right side of equation (7) 
establishes on the left side the validity of the Bernoulli equation for 
the gas flow. 

For a flow which is free from vortices but variable with time, the 
value of the integral is at every instant, independently of the path, 
equal to the potential difference $2 - ®p. However, since the velocity- 

distribution varies with time, the potential difference also becomes t 
dependent on the time, and one obtains the well-known relationship: 


±0 + 


I N 


2 2 

+ V2 + W2 


= 2 )- 




or 


i + -(u 2 + v 2 + w 2 ^) 


+ 


— = Constant 
dt 



This is the generalized Bernoulli equation for the unsteady gas flow 
free from vortices. 

A further application of the induction law is possible for nonsteady 
gas flows which vary periodically with time. In this case, it is easier 
to determine, instead of the. pressure rise at every instant, the mean 
value of the pressure rise for the period of the duration T. If one 
integrates the right side of equation (7) with respect to t over the 
period and then divides by the duration of this period, one obtains the 
time average: - 



dx + v dy + w 


dz) 


J u o 


+ 



dx + v dy + w dz) 


t 0 +T 


( 9 ) 
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Since the velocity field is the same at the times t Q and t Q + T, the 

mean' pressure rise indicated above represents the circulation in one 
instantaneous field on a closed line from Pp via P 2 back to Pq 

divided by the time T. 

Figure 2 shows how to apply this rule to the periodic flow in a 
rotor. One recognizes that the circulation required by equation ( 9 ) 
matches exactly the circulation T around the blade. In the case of 
m blades and an angular velocity of the impeller ca, the period is giv en 

by T = accordingly, the pressure rise of the rotor amounts to: 

mm 


h = L = DSffi ( 9 a) 

T 2jt 

Thus, the significance of the circulation around the blades of the 
impeller with respect to the pressure rise is established for al 1 gas 
flows without increase in entropy. However, it remains to be investi- 
gated whether also the torque and therewith the power consumption 
depends in the same manner on the blade circulation. 


II. TORQUE OF THE IMPELLER WITHOUT MASS FLOW 


The torque of the Impeller in plane flow may be determined from the 
difference of the moments of momentum over a control circle outside the 
impeller and one inside the impeller. If one forms in polar coordinates 
the velocity components in the direction of the radius w r and of the 

circumference w u , with w u counted positive in the direction of the 

increasing angle i|r, there results on the arc element r df of a circle 
of the radius r the mass flow pw r r dty. The momentum of this q uant ity 

in circumferential direction is obtained by multiplication by w u , the 
moment of momentum by multiplication by rw u . By integration of the 
moment of momentum over the entire circumference one obtains a torque 

D = t 2 J pw p w u at (10) 

If the entering fluid in the interior of the impeller is supplied 
without rotation, a sufficiently small radius may exist on which circum- 
ferential components of the velocity do not yet appear. In this case, 
the integral over the external circle according to equation ( 10 ) already 
yields the torque of the impeller. For incompressible flow with 
p = const., the velocity components, w r and w u do not correlate in 
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the absence of obstacles outside of the impeller, so that for the integra- 
tion over v u the mean value of w r may be taken out from under the 

integral sign. The mean value of w r multiplied by the density p Mid 

the circumference of the circle 2rjt, however, is simply the mass G dis- 
charged per second through the impeller. Therefore, in the case of incom- 
pressible flow , the torque may be split into the following factors : 

d = — f i 2 ” v m - <ak = (01) 

2jt Jo Sit ay 

Here m is the number of blades, and F the circulation around each 
of the m blades; the latter transformation results from equation (9a). 
Thus, for an incompressible flow a discharge as well as a circulation 
around the blades is necessary when a torque is to occur. 

Since', according to the results of the previous section, the pres- 
sure rise also depends directly on the blade circulation for compressible 
flow, there would also be a certain justification for assuming that in 
the case of compressible flow no correlation between discharge and blade 
circulation enters equation (10). However, erne single example where a 
torque occurs without discharge or without blade circulation will be 
sufficient for deciding this question in the negative sense. In the _ 
case of the example selected, there appears neither a discharge nor a 
blade circulation, and yet one obtains at higher velocities a torque 
different from • zero . 

In order to represent the flow for the case of vanishing discharge 
and vanishing blade circulation at a large distance ^f ram the impeller, 
one can replace the impeller by a rotating wavy cylinder. The simplest 
wavy cylinder has a cross section in which on the circumference of a 
circle of the radius R, m sinusoidal waves with the wave amplitude A 
are superimposed (fig. 3 for the case m = 3)* In case one wants to 
represent the effect of an impeller with m blades more accurately, one 
could still add further waves with the amplitudes Aq, Aj, etc., and 

the numbers 2m, 3ai, etc. If the amplitude A, for the sake of further 
simplification of the calculation, is limited to small values compared 
to the radius R or, more specifically, compared to the wave 
or 

length L = , one can neglect in equation (8) the square of the gas 

m — . — 

velocity compared to the other s umman ds - Considering constant entropy 
according to equation ( 2 ) one then obtains 


i + M= = Constant 

St J p St 


( 12 ) 
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From the potential <t(r,i|r) "the velocity components are obtained 
in the following manner: 


w T 


_ 

Sr 


w. 


u 


= ia® 

r St 


(13) 


Expressed by these components the continuity equation of the unsteady 
gas flow reads: 


Sp + P 
St r 


|d(rw r ) t ^uj 

dp 

— i * 

d 2 ^ t 1 d3> 

_1_ d 2 ^ 

|_ dr dtj 

- — T P 

dt 

Sr 2 r Sr 

r 2 


= 0 


(Ifc) 


From equations (12) and (14-) there results the well-known equation of 
sound propagation: 


AO = ^ + lM+i S^£ = 1 _ c££ 
Sr 2 r Sr r 2 ^2 a 2 St 2 


(15) 


in which a Is the sonic velocity of the gas according to the following 
relationship : 

a 2 = ^ ( 16 ) 

dp 


Since, for small velocities, the pressure and the density deviate only 
little from the values of the gas at rest p Q and p Q , a may be 
regarded as constant and equal to the value of the sonic velocity for 
this state. By integration of the pressure at constant density, there 
results from equation (12) the pressure 


P = Po ’ P °ft (1?) 

The impenetrability of the surface of the wavy cylinder furnishes the 
boundary condition for the differential equation (15)* At the time t = 0 
the cylinder has the following radii r depending on the central angle ijr 


r = R - A sin rmjf 


( 18 ) 
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Due to the rotation at the angular velocity cu (fig. 3 )> one obtains 
the .dependency of the cylinder radii on i|r and t 


r(i|r,t) = R - A sin m(^ - cot) 


( 19 ) 


Hence, the spinning wavy cylinder produces a radial pumping motion with 
the velocity 


w„ 


= = 

8t 


Pmx> cos m(t - act) 


Because of the smallness of the amplitude A, it is sufficient to pre-. 
scribe this value for the radial component of the gas velocity. According 
to equation (13) , there results the boundary condition for 


w 


r 



= Ahwj cos m(\|r - tut) 


( 20 ) 


Since the integration proper of the differential equation for the 
propagation of sound is known and we are here concerned only with the 
application to the impeller, the detailed calculation and the determina- 
tion of the integration constants may be omitted after stating that the 
elimination of the constants was performed in such a manner that in the 
asymptotic development of the solution for large radii only outgoing 
waves (no incoming waves) were retained. The solution which thereby 
became unique may be written in terms of the Bessel functions of the 
first and second kind J m and Y m and reads _ ■ 
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The phase angle 6 in equation (21), though in itself unessential, 
has the value 


v(=f) 

tan S = i (22) 

m V a I 

If one substitutes this solution into the equation (10) for the torque, 
its integration can be achieved with the aid of the formula-*-: 


Jm( x ) Y m'( x ) - YmteJJip’Cx) « — 

3 tx 


One obtains thus the torque 


2p Q ma 2 A 2 



CVI 

v , /Wr\ 

2 

y [— jj 

+ 

[v jj 



(23) 


The power consumption E = Efco corresponds to the radiated sound output. 


In order to represent the results in dimensionless form, we shall 
use as the Mach number M of the gas flow the ratio of the circumferen- 
tial velocity u and the sonic velocity a: 


M = — = — (24) 

a a 

A coefficient for the resistance to the motion c^ equal to the 

"L P 

torque D divided by the dynamic pressure q = ip Q u , the generating 
surface of the cylinder F = 2Rit, and the radius R is introduced: 


°w “ 


D 


1 

5 * 


,u 2 2R 2 jt 


(23) 


^-Compare Frank-Mises: Dlfferentialgleichungen der Physik (Differ- 

ential equations in physics), 1930 , vol. 1, p. 4l4. 
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and. a second coefficient is formed for the radiated energy c e equal 
to the radiated energy E = Dm divided hy q, F, and the sonic 
velocity ai: 

c e = c^M ( 26 ) 

—p 0 ura2Bjt 


In this form the results of the calculation are plotted in figures 4 
and 5 with the Mach number used as the abscissa and the coefficients Cy- 

or c e used as. the ordinate . The values m indicated at the individual 

curves signify the number of waves on the circumference of the circle. 

The ordinate has as its unit a value which is formed from the amplitude A 
and the wave length L, or the amplitude A, the number of waves m, and 
the radius R. 

One recognizes from the figures that the torque and likewise the 
radiated energy disappears only for incompressible flow and for gas flow 
with very small velocities . In the subsonic domain there results for 
growing m an ever increasing region in which no noteworthy torques 
occur. The maximum of c^. always lies at the Mach number 1. 

It is interesting to compare the rotating wavy cylinder in a resting 
gas and the resting wavy cylinder with a circulatory flow treated by 
G. J. Taylor 2 . While it was found there that even at velocities higher 
than sonic velocity damped perturbation waves occur, the above calcula- 
tion, inversely, yields the result that even below sonic velocity the 
disturbance extends to infinity. Both cases agree at m = °° (that is, 
for a cylinder radius large compared to the wave length) with the solu- 
tion for a flat plate with sinusoidal waves treated before by J.’ Ackeret?. 


SUMMARY 


It is shown by the example of the plane impeller that in a gas 
flow with constant entropy as in an incompressible flow, the pressure 
head of a rotating impeller depends only on the circulation around the 
blades. In contrast to incompressible flow, however, one obtains for 
the impeller rotating freely in an infinite gas mass a larger torque 
than would be necessary for production of the pressure rise because, 
due to the periodic disturbance caused by the rotating impeller, sound 

2 G. IC J. Taylor: ZAMM 10, 1930, p. 334. 

3j. Ackeret: Helvetia Phys. Acta 1, 301j 1928- 
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waves of f ini te energy travel away to infinity. This energy produces a 
resistance which grows with a high initial power of the circumferential 
velocities and, when sonic velocity is exceeded, gradually "becomes the 

3 

wave resistance of bodies moved rectilinear ly at supersonic velocity. 


Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 


See footnote 3 on p. 12. 








